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Abstract: Deep displacement observation is one basic means of landslide dynamic study and 
early warning monitoring and a key part of engineering geological investigation. In our 
previous work, we proposed a novel electromagnetic induction-based deep displacement 
sensor (I-type) to predict deep horizontal displacement and a theoretical model called 
equation-based equivalent loop approach (EELA) to describe its sensing characters. However 
in many landslide and related geological engineering cases, both horizontal displacement and 
vertical displacement vary apparently and dynamically so both may require monitoring. In 
this study, a II-type deep displacement sensor is designed by revising our I-type sensor to 
simultaneously monitor the deep horizontal displacement and vertical displacement 
variations at different depths within a sliding mass. Meanwhile, a new theoretical modeling 
called the numerical integration-based equivalent loop approach (NIELA) has been proposed 
to quantitatively depict II-type sensors' mutual inductance properties with respect to 
predicted horizontal displacements and vertical displacements. After detailed examinations 
and comparative studies between measured mutual inductance voltage, NIELA-based 
mutual inductance and EELA-based mutual inductance, NIELA has verified to be an 
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effective and quite accurate analytic model for characterization of II-type sensors. The 
NIELA model is widely applicable for II-type sensors' monitoring on all kinds of landslides 
and other related geohazards with satisfactory estimation accuracy and calculation efficiency. 

Keywords: electromagnetic induction-based deep displacement sensor; theoretical 
modeling; deep horizontal displacement; deep vertical displacement; mutual inductance 



1. Introduction 

Landslides occur in many areas in the World, causing not only heavy property losses but also 
serious loss of human lives. Landslide deformation is an integrated reflection of geological structure 
and internal and external influencing factors of a landslide mass. Therefore, landslide deformation 
monitoring is a basis for the analysis of the geological structure and dynamic deformation characteristics 
of the investigated landslide mass, a support for informatization design of sliding remediation projects, 
and a promise of feasible technology to predict and provide advance warning against geo-hazards [1-5]. 
Landslide deformation monitoring is based on displacement measurement information and mainly 
includes surface displacement monitoring and deep displacement monitoring [6-8]. Compared to the 
former, landslide deep displacement monitoring is more complicated while more significant, because 
through it, the landslide deformation mode could be deduced, the sliding plane location and depth 
decided, and the dynamic ranges and trends of sliding deformation judged, thereby providing a 
scientific basis and reliable information for the analysis of a landslide's stability conditions, deformation 
mechanics and the related design of treatment engineering [9,10]. 

On the other hand, both surface displacement monitoring and deep (or subsurface) displacement 
monitoring have two aspects: horizontal displacement monitoring and vertical displacement monitoring. 
The significance of measuring and monitoring these two aspects simultaneously has been fully 
demonstrated and verified by lots of theoretical research and engineering cases concerning landslides 
and other slope movements. 

Surface displacement monitoring [11-13], whether in terms of surveying techniques or monitoring 
instrumentation has developed rapidly and reached a high level. The conventional instrumentation for 
surface displacement monitoring include settlement gauges, precision levelings, theodolites, surface 
extensometers for surface vertical displacement measurement, ultrasonic or laser distance meters, and 
deflectometers for horizontal displacement measurement, and total station, aerial photogrammetry for 
measurement in both directions. The modern techniques [14-16] for surface displacement monitoring 
include multi-antenna GPS receiver, DInSAR (Differential Interferometric Synthetic Aperture Radar), 
terrestrial laser scanning, etc. 

Compared to the surface displacement monitoring, development of deep displacement 
monitoring [17,18], whether in terms of monitoring techniques, methods or instrumentations (including 
related sensors) is much slower thus has greatly limited its application scope and popularity, due to 
the extremely complicated and variable characteristics of deep rock and soil mass, such as 
temporal-paroxysmic, spatial-randomicity and invisibility, conditional-terribleness (e.g., corrosion, seepage 
of underground water, failures of shear and compression), geological-heterogeneity and complexity. 
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At present, methods/instruments for deep displacement monitoring can mainly grouped to three 
categories. The first one is borehole extensometers, which includes two subcategories: multi-point 
borehole extensometers [19,20] and wire/cable extensometers [21]. Although wire/cable extensometers 
are relatively simple and low-cost devices, displacements measured by them are global (i.e., the total 
changing distance from one point on landslide surface to another fixed point inside the borehole below 
the slip surface) and the device can neither detect the vertical components and horizontal components 
of the underground slope movement, nor identify the presence of several slip surfaces, nor obtain the 
relative displacements at different depths within the sliding mass. Multi-point borehole extensometers 
are conventional devices to monitor the change of vertical displacements at pre-selected depth that 
decided by pre-installed targets along the common axis of vertical borehole, so they are mainly used 
for settlement and heave monitoring of underground soil and rock. The main disadvantages of borehole 
extensometers include the fact that instrument installation is difficult and complex under deep borehole 
conditions, sliding surfaces are hard to determine, rods/probes may be wedged if large lateral 
displacement occurs, and data reading is laborious. 

The second one is slope inclinometers [22-24]. Presently, they are widely applied to monitor deep 
horizontal displacement at a constant interval of depth within the slope mass and to locate the potential 
sliding surfaces. Although slope inclinometers work based on a relatively simple principle, they do not 
allow continuous recording of displacement, so it is hard to monitor landslides in real-time and 
automatically. Furthermore, they are vulnerable to fail or "shear off due to jams, S-shape or over-large 
displacements, cut-off of inclinotubes and other reasons. 

The third one is Time Domain Reflectometry (TDR) [25-27], which uses an coaxial cable as a 
sensor and works like radar to locate the depth of shear planes or deformation zones in a landslide. 
TDR is a relatively low cost and time-saving monitoring approach, but it can hardly decide the sliding 
magnitude and direction of deep displacement, nor monitor landslides with heavy sliding bands. 
Furthermore, TDR cannot be used where a shear zone does not occur but monitoring tilt is necessary. 

To summarize, at present there exist few sensors or instruments that can simultaneously and 
efficiently monitor the horizontal displacement and vertical displacement of subsurface deformation at 
different depths within the soil and rock mass. 

In our preceding study [28], a novel electromagnetic induction-based integrated deep displacement 
sensor (hereafter called the I-type deep displacement sensor) was presented. It can directly convert the 
varied deep horizontal displacement and tilt angle at any depth within the sliding mass to the 
corresponding variation of mutual inductance between any two adjacent solenoids (each solenoid 
functions as a sensor unit), so it can monitor the underground horizontal displacement more 
conveniently and accurately with a relatively simple and low-cost design. 

Theoretical modeling is an essential and important work in sensor design, error analysis, and 
optimization processes, as it can greatly help to understand the behavior of the sensor so as to optimize 
sensor design and solve some concurrent problems. 

In our previous work [28], in order to describe the complicated relationship among the underground 
sliding mass' horizontal displacement and tilt direction, the I-type sensor's geometric parameters and 
its related output of mutual inductance and mutual inductance voltage, we have initially proposed a 
theoretical model called equation-based equivalent loop approach (EEL A). Through a series of 
comparative studies between the experimental results based on our I-type sensor prototype and 
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EELA-based theoretical simulation results, we not only initially showed the sensor's feasibility, but 
also validated that the EEL A model is quite suitable for depicting the said sensor's sensing properties 
thanks to its commonality, effectiveness, and adequate accuracy. 

We note that in our last paper [28], to adapt to the I-type sensor, it was assumed that along with the 
deep sliding movement, any two adjacent sensor units (labeled as Solenoid I and II) were free to 
relatively tilt and move horizontally, but no obvious vertical displacement occurred between them. The 
same assumption has been adopted in the EELA modeling process for an I-type deep displacement 
sensor. However the supposition that vertical displacement does not change or changes negligibly places 
great limitations on the deep displacement monitoring process of landslides and other geological 
disasters. That's to say, in many monitoring cases for landslides and other slope movements, both the 
horizontal displacement and vertical displacement vary apparently and dynamically, and may require 
simultaneous monitoring of both, which can then provide more comprehensive and objective 
monitoring guidelines for the deep displacement monitoring process. 

Let's explain this in more detail. From the definition of landslide [29-31], both a generalized and 
narrow sense of landslide exists. The general sensor of "landslide" refers to "the movement of a mass 
of rock, debris, or earth down a slope" and mainly includes five types of mass movements: fall, topple, 
spread, slide and flow. The narrow sense refers only to "slides", which mainly include two types: 
translational slides and rotational slides. Translational slides are mainly planarly displaced along the 
sliding surface, so it is reasonable to assume that no obvious vertical displacement occurs inside the 
sliding mass and to mainly monitor the horizontal displacement variations during the deep movement 
process. However, for a rotational landslide, especially during the compression creep stage, both 
vertical displacement and horizontal displacement change obviously and dynamically, so both may 
require monitoring. Moreover, for other landslide types and some related geohazards, especially for 
those caused by excessive underground coal mining, excessive groundwater extraction and slope foot 
cutting, a large amount of theoretical studies and engineering practices show that a simultaneous 
monitoring of deep horizontal displacement and vertical displacement is often required. 

Therefore in this paper, some effective revisions have been made to the I-type deep displacement 
sensor both in the structure design and the correlated driver software, to make it meet the need for 
simultaneously monitoring the horizontal component and vertical component of landslide deep 
displacement. The revised version is then called a II-type deep displacement sensor. 

In order to depict the mutual inductance properties of the proposed II-type sensor efficiently and 
accurately, a theoretical modeling called numerical integration-based equivalent loop approach 
(NIELA) is presented. Combining the numerical integration technique with the equivalent loop 
approach, this model can qualitatively and accurately evaluate the complicated relationships between 
the mutual inductance, the geometrical parameters of any two adjacent sensor units, and their relative 
position in terms of relative horizontal displacement, vertical displacement and axial angle, thereby 
both the variations of deep horizontal displacement and deep vertical displacement together with tilt 
directions at various depths within the monitored slope mass can be simultaneously monitored and 
quantitatively determined by a II-type deep displacement sensor. Modeling verification through 
experimental tests and comparative studies have confirmed the proposed NIELA model's theoretical 
reliability and estimation accuracy in depicting the proposed sensor's sensing properties. 



Sensors 2012, 12 



237 



2. Architecture and Principles of II-Type Deep Displacement Sensor 

As previously described [28], an I-type deep displacement sensor is mainly constituted of a number 
of deep displacement measuring integrated sensor units in series (hereinafter referred to as sensor unit). 
Each sensor unit has identical structure, including an air-core solenoid as main component and 
embedded along its inner wall an integrated sensing circuitry PCB with such functions as sinusoidal 
voltage generation (£//), mutual inductance voltage measuring (U 0 ) 9 tilt angle measurement (#o), RS485 
bus communication with the deep displacement measuring central processing unit, etc. 

Each sensor unit is vertically spaced a certain distance and encapsulated in a heat-shrinkable plastic 
soft tube, so forming a deep displacement measuring chain. By measuring the relative horizontal 
displacement and tilt angle variations between any two adjacent sensor units one by one, the 
cumulative deep horizontal displacement and sliding direction, from surface to different depths within 
the monitored sliding mass can be measured. Any two adjacent sensor units (Solenoid I and II) 
constitute a relatively deep displacement sensor that can be used for measuring the relative horizontal 
displacement and sliding angle at some given depth within the sliding mass. 

Before an I-type sensor takes effect, these sensor units should be vertically buried into a borehole 
and backfilled tightly with cement grout so it can deform synchronously with the surrounding soil 
mass. It is worth noting that to fit the sensing properties of the I-type sensor [28], it is assumed that 
along with the deep sliding movement, any two adjacent sensor units are free to relatively tilt and 
move horizontally, but no obvious vertical displacement occurs between them. From the laws of 
electromagnetic mutual induction, when a sine voltage Ut with fixed frequency and amplitude applied 
to solenoid I, a corresponding mutual induced voltage U 0 and mutual induction M will be generated 
across solenoid II (M has proved to be directly proportional to U 0 ). Under the above assumptions, U 0 
and M have a definitive functional relationship with the relative horizontal displacement X and axial 
angle 9o between solenoid I and II, and their geometrical parameters in terms of diameter D, length A 
and winding coil turns W. So, the functional relation can be generally expressed as: 

U o =f{X,0„Z„D,Ajy) (1) 

where Z 0 is the initial axial distance between Solenoid I and II, Oo is their relative axial angle (equaling 
to the relative tilt angle) which is measurable by the sensor' s integrated tilt measuring unit. 

In order to describe the above mentioned complicated relationship [Equation (1)] analytically and 
accurately, in our former study [28], we have initially proposed the EELA model for I-type sensors. 
EEL A has been tested to be a reliable and effective model to depict I-type sensors' sensing properties. 

However just as previously explained, it is a considerable limitations to suppose the vertical 
displacement does not vary or varies negligibly in the sliding deep displacement monitoring process, 
whether in theory or on practice. Therefore in this study, a II-type deep displacement sensor is 
proposed by making revisions to the structural design and the corresponding driver software of the 
I-type sensor. 

Briefly speaking, the main structural revisions to each sensor unit of a II-type sensor includes a 
small cylindrical permanent magnet mounted at the center of the lower surface and a high sensitivity 
linear Hall sensor located at the center of upper surface (Figure 1). According to the Hall effect, the 
output voltage of Hall sensor in magnetic field satisfies the following change rule: 
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U H =R H I -^- = R H B (2 ) 

where R H and d are the Hall coefficient and semiconductor slice thickness of the Hall sensor 
respectively, In is the current applied on the Hall sensor, B is the magnetic induction intensity exerted 
in the direction perpendicular to the upper surface of the sensor package. 

Figure 1. Structural diagram of a II-type deep displacement sensor unit. 




As Figure 2 shows, when Solenoid I and II produce a relative displacement (whether horizontal 
displacement or vertical displacement or a combination of both), the relative position changes between 
the permanent magnet on the upper surface of Solenoid I and the Hall sensor on the lower surface of 
Solenoid II, so the magnetic field applied on the Hall sensor is changed accordingly, thus making the 
correlated Hall sensor output voltage change. That is, there exists a certain functional relationship 
between the magnetic induction intensity B and the relative position of Solenoid I and II. Considering 
the axial symmetry of cylindrical permanent magnet, the generated magnetic field also shows axial 
symmetry, which may be described as follows: 

B = f l (X,Z,0 o ) (3) 

where X and Z are the central distance and axial distance between Solenoid I and II and can be 
used to describe the relative horizontal displacement AX and relative vertical displacement AZ, 
respectively, i.e., 

AX = X-X Q 

AZ = Z-Z 0 ^ 

where Xq and Zo are the initial central distance and axial distance respectively between these two 
solenoids, generally setting Xq = 0 and Z 0 = 40 mm. 
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Figure 2. Geometrical arrangement between two arbitrary adjacent solenoids: (a) Initial 
arrangement; (b) Relative horizontal displacement and vertical displacement occurred; 
(c) Relative tilt, horizontal displacement and vertical displacement occurred. 




(a) (b) (c) 

There are some methods [32-35] to calculate the magnetic field B of a cylindrical permanent 
magnet, including the equivalent magnetic charge method, equivalent electric dipole method, 
equivalent current method and finite element method. Here we wouldn't detail these due to the paper 
length limitations. 

Combining Equations (2) and (3), the relationship between Solenoid IPs output Hall voltage Uh and 
its position relative to Solenoid I can be described as: 

U H =tiMX,Z,0 o ) (5) 

Meanwhile, the functional relationship among Solenoid IPs output mutual inductance voltage U 09 
the relative geometrical position between Solenoid I and II in terms of relative horizontal displacement 
(AX = X - Xo), vertical displacement (AZ = Z - Zo) and axial angle (#o), and their geometrical 
parameters in terms of diameter D, length A and coil turn Wean be expressed as: 

U o =f 2 (X,Z,0 o ,D,A,W) (6) 

Combing Equations (5) and (6), it shows that a II-type sensor requires no assumptions of no relative 
vertical displacement taking place between two adjacent solenoids during deep sliding monitoring 
process, so for any two adjacent sensor units, whether variations of relative horizontal displacement, 
vertical displacement, or axial angle, can cause the mutual inductance voltage U Q and the Hall sensor 
output voltage Unto change at the same time. During the working process of II-type sensor, so long as 
U 09 Uh and 9o between any two adjacent sensor units can be synchronously and automatically 
measured by the proposed sensor itself, and Equations (5) and (6) can be accurately expressed and 
evaluated by theoretical modeling or equation derivation, the said sensor can quantitatively determine 
the relative horizontal displacements and vertical displacements at different depths within the sliding 
mass. It is worth mentioning that the probable measured results include two special conditions: (i) the 
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landslide is totally caused by horizontal displacement, then the sensor will measure AZ ^ 0; (ii) the 
landslide is totally caused by vertical displacement, then the sensor will measure AX ^ 0. 

As mentioned above, the evaluation of Equation (5) is relatively simple with some available models 
and solving methods for reference. Compared to that, the evaluation of Equation (6) is much more 
complicated and remains a difficult issue for which there are few existing evaluation equations or 
models nowadays, so in this paper, one of the main theoretical tasks is to develop an efficient and 
accurate theoretical model to depict the complex functional relationship among mutual inductance 
voltage U 09 the geometrical parameters of Solenoid I and II, and their relative central displacement X, 
axial distance Z and axial angle #o, which reflect directly the relative sliding horizontal displacement, 
vertical displacement and tilt angle at the sliding mass's corresponding depth. 

It is worth stressing, just as our previous work has shown, that although both the change of mutual 
inductance voltage U Q and that of mutual inductance M respond to the variations of relative 
displacement and axial angle between Solenoid I and II, it's much simpler and more efficient to 
investigate deep displacement in terms of mutual inductance rather than mutual inductance voltage. 
Meanwhile, M is strictly proportional to U 0 and their relationship can be expressed as: 

TT 1 dU i M 

U= -M n\ 

0 R dt K } 

where R is the equivalent resistance of Solenoid I, and Ut is the 10 KHz sine voltage applied on 
Solenoid I. 

Therefore, in this paper, we focus on II-type deep displacement sensor, and use the general 
Equation (8) to depict the functional relationship among mutual inductance M, the geometrical 
parameters of Solenoid I and II in terms of diameter D, length A and winding coil turns W, and their 
relative position in terms of central distance X 9 axial distance Z and axial angle #o: 

M = f 3 (X,Z,0 o ,D 9 A,W) (8) 

In order to evaluate the above complicated relationship [Equation (8)] qualitatively and effectively, 
we present in this paper a theoretical model called NIELA after extensive and intensive researches. 

Compared to the existing EELA model, the proposed NIELA model uses a numerical integration 
approach rather than an equational derivation in infinite series form to evaluate the mutual inductance 
M , to meet the modeling requirement of varying both relative vertical displacement and horizontal 
displacement. Therefore NIELA is applicable not only for II-type sensors but also I-type sensors, 
whereas, EELA, it is generally only applicable for I-Type sensors, because the infinite series 
expressions for mutual inductance may be unconvergent and become invalid when varying the relative 
vertical displacement between Solenoid I and II. 

3. NIELA for II-Type Deep Displacement Sensor 

EELA was previously introduced in detail [28]. Here we only present a summary of this approach 
to explain why EELA is suitable for I-type sensors but not for II-Type sensors. Then we will introduce 
in detail NIELA, which is applicable both for I-type and II-type deep displacement sensors. 

In brief, the EELA model is essentially an approximate calculation based on the double integration 
with the following basic steps: 
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Step 1: using two equivalent current loops to replace one solenoid, so the mutual inductance 
between any two adjacent solenoids can be equivalent as: 

M = W 2 (M 13 + M 14 + M 23 + M 24 ) / 4 (9) 

where M 13 , M 14 , M 23 and M 24 are the mutual inductances between two equivalent loops 1 and 3, 1 
and 4, 2 and 3, 2 and 4, respectively. 

Step 2: applying some related electromagnetic field theory and equations to deduce the equational 
expressions of mutual inductance for M 13 , M 14 , M 23 and M 24 , respectively. 

Figure 3. Geometric arrangement between equivalent parallel-axial Loop 2 and Loop 3. 
(a) Front view, (b) Top view. 
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For example, when Solenoid I and II are in parallel-axial arrangement, then M 13 , M 14 , M 23 and 
M 24 are the mutual inductances between two equivalent parallel-axial current loops 1 and 3, 1 and 4, 2 
and 3, 2 and 4, respectively. When Loop i and Loop j are arranged as Figure 3 shown, we can deduce 
the following mutual inductance expression for M .. : 



m =M5yr-ir 1 — 
9 i6 tr n+l 



(2rc-l)!! 
(2n)!! 



So long as the following convergence condition be satisfied: 

1, < 1 



(10) 



(11) 



where i = {1,2} , y = {3,4} , R=R=D/2, r^^Xl+Zl , A = Z>/r,, fy=Z^, P 2w (/7,) is the 
Legendre polynomials of degree 2n with argument 77. , and jU 0 =47rxl0~ 7 H / rn is the free 
space permeability. 

Obviously, the expression for M fj is quite complicated and expressed by an infinite series, which 
means, when X tj « 1 , M.. is dominated by the first several terms because the series converges quickly, 
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but the more Ay is close to 1, the more slowly the series converges, and the more terms need be 
calculated to get an approximation to M.. with sufficient accuracy, so calculation of M 14 is very 
time-consuming when Ay is very near 1. Much worse, when Ay happens to be 1 or larger than 1, then 
Equation (10) cannot converge at all. 

From Equation (1 1), it is seen that, under a given set of conditions, A ij will be smaller if r tj is larger. 
When Solenoid I and II are in parallel state, X y is equal to X, which is determined by the relative 
horizontal displacement between Solenoid I and II that occurred along with the slope movement. 
That's to say, when the horizontal movement of slope is very slow or the sensor is buried in the slope 
mass not long ago, X y may be so small (for instance, Xy = 5-10 mm) that r tj is mainly determined by 
Z tj . Note that in these four equivalent loops, Loop 2 and 3 are the closest to each other, so Z 23 is the 
smallest one of Z^.(Z 13 ,Z 14 ,Z 23 ,Z 24 ), that is, if Z 23 satisfies the convergence condition, then the other 
three Z.. are sure to converge. Now we will examine Z 23 in detail. 

For the experimental sensor prototype, the solenoids' diameter and length are set to be D = 70 mm 
and A = 75 mm. According to the equivalent loop approach, approximately Z 23 - Z-AI 4^ . Then the 
convergence condition for Z can be expressed by: 



For example, if X = 5 mm, Z must be larger than 113.1 mm; if X = 10 mm, Z must be larger 
than 1 12.6 mm, and so on. 

So for an I-type sensor and the correlated experiments conducted before [28], the initial value of Z 
was set as 1 15 mm and supposed to not vary with the sliding movement. Under such an arrangement, 
and all other Ay could be guaranteed to converge, so Equation (10) could be quickly convergent 
and thus effective in calculation. From this it follows that EEL A is applicable to I-type sensors. 

Furthermore, after a series of comparisons and examinations conducted [28] between the predicted 
mutual inductance based on EELA and the experimentally measured mutual inductance voltage based 
on an I-type sensor prototype, the EELA model was tested to be reliable and effective in depicting an 
I-type sensor's sensing properties (i.e., determining the relative horizontal displacement and tilt angle 
quantitatively) with acceptable estimation accuracy on the premise of convergence. 

However, considering there exist limitations in assuming the vertical displacement does not vary in 
the sliding process whether in theory or on practice, the II-type deep displacement sensor is proposed 
to free us from this assumption. That is, for a II-type sensor, any two adjacent sensor units are free to 
relatively tilt (#o), move horizontally (Z), and move vertically (Z) along with sliding of the surrounding 
rock and soil mass. 

Under such instances, when Z is reduced from the initial 115 mm to 1 10 mm or less, Equation (10) 
will no longer converge and becomes invalid to evaluate M.. . When Solenoid I and II are arranged in 
cross-axial state, the equational expressions deduced for M tj are also complicated and form an infinite 
series, which we have derived in detail before [28]. All these facts show that EELA is basically invalid 
to depict II-type sensors due to the non-convergence problem during variation of relative vertical 
displacement between any two adjacent sensor units. 

In this paper, a new theoretical modeling named numerical integration-based equivalent loop 
approach (NIELA) is proposed to depict the mutual inductance properties of II-type sensor. This model 
can qualitatively depict the complicated relationships among mutual inductance M, geometrical 




(12) 
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parameters of Solenoid I and II in terms of diameter D, length A and coil turns W 9 and their relative 
position in terms of relative horizontal displacement (AX=X-Xo), vertical displacement (AZ = Z - Z 0 ) 
and axial angle (Oo) for any two adjacent sensor units just as Equation (9) denoted. 

Compared to EELA, NIELA applies the same hypotheses for the modeled air-core solenoids [28] 
and the same equation [i.e., Equation (9)] to evaluate the mutual inductance between any two adjacent 
solenoids, but uses numerical integration rather than infinite series to express and evaluate M u , M l4 , 
M 23 and M 24 . For convenience of interpretation, we will demonstrate how M 14 is evaluated in the 
NIELA model when Loop 1 and 4 are arranged in a parallel-axial state as shown in Figure 3. 

Firstly, the Cartesian coordinate O-xyz and polar coordinate O-pcpz are established simultaneously, 
in which we let Loop 1 and Loop 4 lie in the xy plane, having radii R x and R 4 , respectively, and apart 
from each other by an axial distance Z 14 and a central distance X l4 . Supposing Loop 1 carries current 
I x , then under the polar coordinate system, an arbitrary source point Q in Loop 1 can be denoted as 
Q{p l ,(p l ,Z l ) = Q(R l9 cp l9 0) , and according to the Biot-Savart law, the vector potential at an arbitrary 
field point P{p 9 (p 9 z) due to current I x in Loop 1 is: 

Ai(A ^ z)= IHir (13) 

In Equation (13), the integration is along the direction of current flow, the current element l x dt is 
tangent to Loop 1 at source point Q , R is the distance vector from the source dt to the field point P , 
R = \R\ and R = R/ R,jU 0 =4^x10"' HI m is the free space permeability. After simplification: 

A l (p,<p 9 z) = A pl e p +A 9l e 9l ( 14 ) 



A 



_ jUqI^ sin (j)d(j) 



sin (pay 

PX ~ ^ >-«Jp 2 +R?-2pR lC os0 + z 2 (15) 

where (j) = (p-cp x . Using parity of trigonometric functions, we can prove the integral of p component 
of Equation (14), A x =0 9 so the vector potential is azimuthal, i.e., A l (p 9 <p,z) = A l (p,z)e l . Let 
0 = 7r-2j3 , then: 

n jift + Ri) 2 +z 2 -4/^sin 2 p 

Jo L , 2 ■ 2 o 



xjip + Rif+z 2 Jo jl-k 2 sin 2 p 

=mVva/w/(2*) 

and: 

f(k) = (2/k-k)K(k)-2E(k)/k (18) 

where k = ^4pR l I [(p + R l ) 2 + z 2 ] , AT and E are the complete elliptic integrals of the first and second 
kinds respectively with modulus k : 

K(k\-[" 2 d & 
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E(k) = \* j\-k 2 sin 2 fidfi (20) 

If we limit P to be one point in Loop 4 as Figure 3 shown, then: 

P(p,(p,z) = P(p 4 ,<p 4 ,z 4 ) (21) 

A = V X i4 +i? 4 +2X 14 R 4 cos#> 4 (22) 

A 2 + ^ 4 2 - X 2 U R 4 + X l4 cos (p 4 

2* 4 A * 4 A 1 JJ 

Z *= Z M (24) 

According to the electromagnetic induction theory, the mutual inductance between Loop 1 and 
Loop 4, M 41 = <E> 41 / /j , where 0 41 is the magnetic flux through Loop 4 due to current I y in Loop 1. 
And 0 41 can be evaluated by: 

*4i=j> C4 MP)-<U* (25) 

In our case: 



A l C4 



M i4= — $ A <pie<pre<p4 dC 4 



1 f 2* , x (26) 

To evaluate M 14 explicitly, we first calculate M 0 , the mutual inductance between two coaxial 
current loops whose radii are R x and p 4 respectively, and Z 14 apart inz axis: 

1 r 27r ~ ~ 27rp d 

M o=yl MP> z >, l ' e , l P d n = —r-MP*> z iJ (27) 

Then M 14 can be associated with M 0 by: 

M u = C[2ffp 4 A <pl (p 4 ,Z u )/I l ][R 4 cosa/(2xp 4 )] d<p 4 

(28) 

= f * M 0 [R 4 cos a/ (2np 4 )]dq> 4 

J 0 

Let <p = (p 4 , then: 

M u = M 41 = — I — ^ <ty (29) 

cos a = (7? 4 + X 14 cos q>) I (R 4 p 4 ) ( 30 ) 

Combining Equation (29) with Equations (17-20), (22), (27), (28), and (30), M l4 can be easily 
evaluated by numerical integration over the range 0<(p<7t . From this, we can see the convergence 
limitations set upon EELA has been completely overcome by NIELA, so NIELA is applicable for the 
proposed II-type deep displacement sensor. 
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4. Experimental Testing and Model Verification 

4.1. Experimental Setup and Procedure 

To verify the above analysis and to test the performance of the proposed NIELA method on 
evaluating of the mutual inductance versus the horizontal displacement and vertical displacement 
between two adjacent sensor units for a II-type sensor, we conducted a series of experiments and 
comparisons using a sensor prototype and some related devices, which include the sensor's axial, 
horizontal displacement and vertical displacement drive devices, axial angle measurement unit, sensing 
data acquisition, processing, communication and display unit, etc. The sensor prototype mainly 
includes two adjacent integrated deep displacement sensor units (Solenoid I and II) and a deep 
displacement measuring central processing unit. Under the control of the central processing unit, the 
sine input voltage Ut can be automatically generated on Solenoid I, the corresponded mutual inductance 
voltage U 0 across Solenoid II and the axial distance Z, central distance X and axial tilt angle 9o between 
them can be automatically measured and recorded in real time and further transmitted to a remote or 
local comprehensive processing center for detailed process through RS-485 or wireless communication. 
X, Z and 9o can consecutively adjusted by the sensor's axial, horizontal displacement and vertical 
displacement drive devices. The detailed sensor fabrication process and supported devices arrangement 
have been introduced before [28] and a photograph of the experimental setup is shown in Figure 4. 



Model verification process mainly includes two parts: (i) test the modeling effectiveness of EELA 
and NIELA for an I-type sensor (X Variable, Z invariant); (ii) test the modeling effectiveness of 
EELA and NIELA for a II-type sensor (both X and Z Variable). This is conducted mainly by 
comparison among the measured mutual inductance voltage, NIELA-based mutual inductance, 
EELA-based mutual inductance under the same given conditions. It is noted that the change of mutual 
inductance should be completely proportional to mutual inductance voltage in theory [28]. 

4.2. Experiments and Model Validation One (Z not varied) 



Figure 4. Photograph of the sensor prototype based experimental setup. 




To test the modeling effectiveness of NIELA and EELA for an I-type sensor (X variable, Z 
invariant), we first conducted experiments using the following assumptions: under the impact of deep 
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slope movement, both the relative horizontal displacement X and tilt angle 6o are changed between 
Solenoid I and II, but their relative vertical displacement (AZ = Z - Z 0 ) does not vary or 
varies negligibly. 

Figure 5. 3-D graphs of (a) measured mutual inductance voltage; (b) NIELA-based 
mutual inductance; (c) EELA-based mutual inductance versus center distance and axial 
angle between Solenoid I and II (Z = 115 mm). 




(c) 
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Table 1. Geometrical parameters of Solenoid I and II. 



Parameter 



Unit 



Value 



Comment 



Diameter (D) 
Length (A) 
Axial distance (Z) 
Coil turns (W) 



mm 



mm 



mm 



mm 



70 
75 
115 
400 



divided by 3 layers 



Obviously, an I-type sensor is workable in such instances, so in the experiment, we fixed the axial 
distance Z to be 115 mm (i.e., Z = Zo) but gradually varied X (0-1 00 mm, range interval: AX= 2.5 mm) 
and 9o (0-75°, range interval: A9 = 5°), and recorded the corresponding output of the mutual 
inductance voltage between these two solenoids, and finally plotted these measured data into 3-D 
graphs, as shown in Figure 5(a). Meanwhile, we plotted the corresponding 3-D theoretical prediction 
graphs based on NIELA and EEL A respectively, which are shown in Figure 5(b,c). The geometrical 
parameters for the modeled and sensor prototype-based solenoids are listed in Table 1 . A comparison 
of Figure 5(b,c) to Figure 5(a) shows that very good agreements are achieved between the experimental 
data and modeling output wherever Solenoid I and II are in coaxial, parallel-axial or cross-axial states, 
which indicates both the NIELA and EELA models are quite reliable and effective to describe the 
property of an I-type deep displacement sensor. 

To allow further examinations, some 2-D curves were extracted from its 3-D graphs in Figure 5 by 
specifying some axial angle 9o (for instance, 5° and 25°) as Figure 6 shows. A comparison between 
Figure 6(a) and Figure 6(b,c) shows that the experimental data still show high shape similarity to 
modeling results based on both NIELA and EELA, thereby further verifying these two models' 
reliability and high approximation in formulation of an I-type sensor's sensing properties under the 
hypothesis that no relative vertical displacement occurred between any two adjacent sensor units inside 
the sliding mass. 

Figure 6. 2-D graphs of (a) measured mutual inductance voltage; (b) NIELA-based 
mutual inductance; (c) EELA-based mutual inductance versus center distance between 
Solenoid I and II (Z = 115 mm). 




(a) 
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Figure 6. Cont. 




(c) 

4.3. Experiments and Model Validation Two (Z varied) 

As we know, to suppose the vertical displacement does not change when landslides and related 
geo-engineering happens does not fully meet the practical situation of sliding movement, so in this part, 
we will study the influence of both the deep vertical displacement and deep horizontal displacement 
(X & Z variable) on the sliding mass and fully examine the modeling effectiveness of NIELA for a 
II-type sensor under such circumstances. A series of comparative experiments were conducted among 
the measured mutual inductance voltage, the predicted mutual inductance based on NIELA and EELA 
respectively, versus the simultaneous variation of axial distance Z and central distance X under some 
fixed axial angers 9o (Oo can be automatically measured by a II-type sensor ) between Solenoid I and II. 

It can be seen, Figures 7 and 8 plot the 3-D graphs of (a) measured mutual inductance voltage, 
(b) NIELA-based mutual inductance and (c) EELA-based mutual inductance versus center distance X 
and axial distance Z between Solenoid I and II under two different conditions, respectively: 

Condition 1:9 0 = 0°,Z = 101-130 mm 
Condition 2: 9 0 + 0°, Z = 101-130 mm 
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Figure 7. 3-D graphs of (a) measured mutual inductance voltage; (b) NIELA-based 
mutual inductance; (c) EELA-based mutual inductance versus axial distance and center 
distance (0 O = 0°, Z= 101-130 mm). 




120 

Axial distance Z (mm; 



130 80 



Center distance X (mm) 



(c) 
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Figure 8. 3-D graphs of (a) measured mutual inductance voltage; (b) NIELA-based 
mutual inductance; (c) EELA-based mutual inductance versus axial distance and center 
distance (0 O = 15°, Z= 101-130 mm). 




120 

Axial distanca Z (mm) 



130 80 



Center distance X (mm) 



(b) 



EELA-based Mutual indutance (Axial angle 6 n =15° ) 




120 

Axial distance Z (mm; 



130 80 



Center distance X (mm) 



(c) 
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A series of comparative studies show that: 

(i) Under Condition 1, where the convergence conditions for EELA cannot be satisfied, a great 
discrepancy has occurred between the measured mutual inductance voltage and the EELA-predicted 
mutual inductance, so EELA is shown to be invalid due to its divergence, whereas, under the same 
conditions, the NIELA-based mutual inductance shows high consistency to the measured mutual 
inductance voltage, so the NIELA model is tested to be valid and effective under such conditions. 

(ii) Under Condition 2, where EELA satisfies the convergence conditions, both the NIELA-based 
and EELA-based mutual inductances show good tracking of the measured results of mutual inductance 
voltage, so NIELA is still verified to be feasible and effective in modeling a II-type sensor under 
condition 2. EELA also seems effective under this condition. 

(iii) To allow further studies, as shown in Figures 9 and 10, from the 3-D graphs in Figures 7 and 8, 
we have extracted some 2-D curves by fixing the value of X 9 which offers a close-up view of the effect 
of axial distance Z on the measured mutual inductance voltage and predicted mutual inductance based 
on NIELA and EELA, respectively, under some specific central distance X Figures 9 and 10 show the 
simulation results in parallel-axial and cross-axial state, respectively. 

It is noted that, only with the premise of convergence could EELA be correctly apply to theoretical 
modeling for the deep displacement sensor. That is, when Z is smaller than a certain fixed value (i.e., 
convergence limit), the convergence condition for EELA could no longer be satisfied, so the predicted 
mutual inductance is meaningless and invalid, which has been clearly demonstrated by Figure 9(c). 
Meanwhile, the NIELA-based 2-D theoretical curves [Figure 9(b)] are seen to be quite in agreement 
with the experimental one [Figure 9(a)] whereverX= 2.5 mm, 20 mm or 35 mm. 

Figure 9. 2-D graphs of (a) measured mutual inductance voltage; (b) NIELA-based 
mutual inductance; (c) EELA-based mutual inductance versus axial distance Z (Oo = 0°, 
Z= 101-130 mm). 




(a) 
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Figure 9. Cont. 




(c) 

Figure 10. 2-D graphs of (a) measured mutual inductance voltage; (b) NIELA-based 
mutual inductance; (c) EELA-based mutual inductance versus axial distance Z (# 0 = 15°, 
Z = 101-130 mm). 




(a) 
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Figure 10. Cont. 




(c) 

Figure ll(a-c) further plot the normalized curves of the measured mutual inductance voltage, 
NIELA-based mutual inductance, and EELA-based mutual inductance respectively according to 
Figure 9. It can be seen that, under such stringent point to point spatial comparison, the NIELA-based 
mutual inductance still shows good tracking of the measured mutual inductance voltage, which further 
increases our confidence in using the proposed NIELA model to predict both horizontal displacement 
and vertical displacement variations for the II-type sensor. However for EELA, so long as the 
convergence conditions are not met, its normalized curves of mutual inductance show little shape 
similarity to the measured one, so the EELA model becomes invalid and unqualified for II-type sensors 
under such instances. As can be seen from Figure 10, under cross-axial state, both NIELA-based and 
EELA-based mutual inductance match the actual mutual inductance voltage very well. In the same 
way, we have further plot the normalized curves from Figure 10 and labeled them as Figure 12, which 
shows that even under such stringent point to point spatial comparison, these three normalized curves, 
that is, NIELA-based predictions, EELA-based predictions and sensor prototype-based measurement 
results, are almost completely overlapped wherever X— 2.5 mm, 20 mm or 35 mm. 
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Figure 11. Normalized curves for measured mutual inductance voltage, NIELA-based 
mutual inductance and EELA-based mutual inductance when 6o = 0°. (a) X = 2.5 mm; 
(b) X= 20 mm; (c) X= 35 mm. 



Measured mutual indutance voltage vs. Modeled mutual indutance ( 8 n = 0° ) 



□ Measured mutual inductance voltage 
Mutual indutance (EELA-based) 
-e- Mutual indutance (HIELA-based) 



( Central distance X=2.5mrn ) 




110 115 120 

Axial distance Z (mm) 



(a) 

Measured mutual indutance voltage vs. Modeled mutual indutance ( 6 0 =0° ) 



□ Measured mutual inductance voltage 
Mutual indutance (EELA-based) 
-©- Mutual indutance (HIELA-based) 



( Central distance X=20mm ) 




(b) 



Measured mutual indutance voltage vs. Modeled mutual indutance ( B Q =0° ) 



□ Measured mutual inductance voltage 
Mutual indutance (EELA-based) 
-©- Mutual indutance (HIELA-based) 



( Central distance X=35mm ) 




110 115 120 

Axial distance Z (mm) 



(c) 
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Figure 12. Normalized curves for measured mutual inductance voltage, NIELA-based 
mutual inductance and EELA-based mutual inductance when #o = 15°. (a) X — 2.5 mm; 
(b) X= 20 mm; (c) X= 35 mm. 



Measured mutual indutance voltage vs. Modeled mutual indutance ( B 0 =15° ) 



□ Measured mutual inductance voltage 
Mutual indutance (EELA-based) 
-e- Mutual indutance (HIELA-based) 



( Central distance X=2.5mrn ) 




110 115 120 

Axial distance Z (mm) 



(a) 

Measured mutual indutance voltage vs. Modeled mutual indutance ( 9 n =15° ) 



□ Measured mutual inductance voltage 
Mutual indutance (EELA-based) 
-O- Mutual indutance (HIELA-based) 



( Central distance X=20mm ) 




(b) 



Measured mutual indutance voltage vs. Modeled mutual indutance ( 6 Q =15° ) 



□ Measured mutual inductance voltage 
Mutual indutance (EELA-based) 
-©- Mutual indutance (HIELA-based) 



( Central distance X=35mm ) 
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This further validates that NIELA models could quite accurately estimate the variations of both 
horizontal displacement and vertical displacement during sliding movement for the proposed II-type 
deep displacement sensor, whether any two adjacent sensor units relatively slide horizontally, 
vertically or tilt. 

5. Conclusions 

Landslides are one of the most costly catastrophic events in terms of human lives and property 
losses. Deep displacement monitoring is one basic means of dynamic study and early warning 
monitoring of landslides. It is also an important part of engineering geological investigation. Presently, 
there are few single sensors or instruments that can simultaneously and efficiently monitor the deep 
horizontal displacements and vertical displacements from surface to different depths within the monitored 
soil and rock mass on purpose of sliding geohazard monitoring or treatment engineering assessment. 

In our previous work, we have proposed an electromagnetic induction-based deep displacement 
sensor (I-type deep displacement sensor) and a corresponding theoretical model called EELA to 
predict the deep horizontal displacement at different depths within the landslide mass. 

In this study, in order to meet the engineering requirement of monitoring both the horizontal 
displacement and vertical displacement at different depths within the sliding mass, the II-type deep 
displacement sensor is proposed by modifying the I-type sensor. Compared to the I-type sensor, 
whether the variations of relative horizontal displacement, vertical displacement, or axial angle 
between any two adjacent sense units, can cause the mutual inductance voltage U 0 (which is proportional 
to mutual inductance M) and the Hall sensor output voltage Uh to vary simultaneously, so a II-type 
sensor need not make assumptions that no relative vertical displacement occurred inside the slope mass. 
In all, the proposed II-type sensor combines deep horizontal displacement and vertical displacement 
monitoring capability. 

To depict a II-type sensor's mutual inductance properties analytically and quantitatively, a 
theoretical model called numerical integration-based equivalent loop approach (NIELA) is presented. 
Combining numerical integration technique with equivalent loop approach, this model can quite 
accurately evaluate the complicated relationship among the mutual inductance, the geometrical 
parameters of any two adjacent sensor units, and their relative position (i.e., horizontal displacement, 
vertical displacement and tilt angle) just as Equation (8) denoted, through which to predict both deep 
horizontal displacement and vertical displacement variations for a II-type sensor. 

To test the NIELA model's theoretical reliability and estimation accuracy for the proposed II-type 
sensor, a series of comparisons and examinations have been conducted between the measured mutual 
inductance voltage, NIELA-based mutual inductance and EELA-based mutual inductance under 
several application circumstances, from which some main conclusions can be drawn as follows: 

(1) "Experiments and Model Validation One" (where Z is assumed to not vary during slide process 
so the convergent condition for EELA can be guaranteed) shows that very good agreements have been 
achieved among the experimentally measured data, NIELA-based predictions and EELA-based 
predictions, which indicates both NIELA and EELA can effectively and quantitatively express the 
sensing properties of an I-type deep displacement sensor. 
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(2) Through "Experiments and Model Validation Two" (where both X and Z varied), we can see 
that: (i) a great discrepancy has occurred between the measured mutual inductance voltage and 
EELA-based predicted mutual inductance due to a lack of convergence for EELA when varying the 
relative vertical displacement between any two adjacent senor units, so the EELA model is tested to be 
basically invalid for the II-type sensor due to the convergence limitations. Secondly, the NIELA-based 
mutual inductance is found to be in good agreement with the measured mutual inductance voltage, 
which indicates that NIELA is a relatively accurate and efficient model to predict both the deep 
horizontal displacement and vertical displacement for the proposed II-type sensor. 

(3) In sum, the NIELA-based predicted mutual inductance always shows good tracking of the 
measured mutual inductance voltage under all conditions in any experiments conducted, even 
including the most stringent point to point spatial comparisons between them. It can be said that all 
experiments conducted here have verified the NIELA model's high theoretical reliability and 
prediction accuracy in depicting of the mutual inductance characters of II-type deep displacement 
sensor, so both the deep horizontal displacement and vertical displacement at different depths within 
the slope mass could be quantitatively predicted. 

These conclusions, in turns, support these two opinions: 

(1) EELA is well qualified to describe the sensing characters of an I-type deep displacement sensor, 
which is mainly applied to monitor such landslides and related geo-engineering whose main form of 
movement is horizontal displacement while the vertical movement is relatively small or unimportant. 

(2) NIELA is a quite reliable and high approximation model to describe the sensing properties both 
for I-type and II-type deep displacement sensors, so it is generally applicable for monitoring of 
different kinds of landslides and some related geo-engineering problems, especially for such 
monitoring circumstances that both the underground vertical displacement and horizontal displacement 
change dynamically during the sliding process thus a simultaneous monitor toward both displacements 
may really required. 

Acknowledgments 

This work is funded by the National Natural Science Foundation of China (NSFC) of the Special 
Fund for Basic Research on Scientific Instruments under Grant 61027005 and the NSFC General 
Project under Grant 51074146. 

References 

1. Moss, J.L.; McGuire, W.J.; Page, D. Ground deformation monitoring of a potential landslide at 
La Palma, Canary Islands. J. Volcanol Geotherm. Res. 1999, 94, 251-265. 

2. Dai, F.C.; Lee, C.F.; Ngai, Y.Y. Landslide risk assessment and management: An overview. Eng. 
Geol. 2002, 64, 65-87. 

3. Kilburn, C.R.J.; Pasuto, A. Major risk from rapid, large-volume landslides in Europe (EU Project 
RUNOUT). Geomorphology 2003, 54, 3-9. 



Sensors 2012, 12 



258 



4. Qi, S.; Yan, F.; Wang, S.; Xu, R. Characteristics, mechanism and development tendency of 
deformation of Maoping landslide after commission of Geheyan reservoir on the Qingjiang River, 
Hubei Province, China. Eng. Geol. 2006, 86, 37-51. 

5. Barla, G.; Antolini, F.; Barla, M.; Mensi, E.; Piovano, G. Monitoring of the Beauregard landslide 
(Aosta Valley, Italy) using advanced and conventional techniques. Eng. Geol. 2010, 116, 218-235. 

6. Corominas, J.; Moya, J.; Ledesma, A.; Lloret, A.; Gili, J. Prediction of ground displacements and 
velocities from groundwater level changes at the Vallcebre landslide (Eastern Pyrenees, Spain). 
Landslides 2005, 2, 83-96. 

7. Wang, F.; Okuno, T.; Matsumoto, T. Deformation characteristics and influential factors for the 
giant Jinnosuke-dani landslide in the Haku-san Mountain area, Japan. Landslides 2007, 4, 19-31. 

8. Lu, M.-C; Tang,T.-Y.; Tsai, C.-P.; Cho, Y.-A. The non-contact landslides monitoring system for 
long-distance. In Proceedings of the 2010 International Conference on Internet Technology and 
Applications, 20-22 August 2010; pp. 1-4. 

9. Jin, X.; Li, X.; Wang, L.; Wang, X. Characteristics of landslide deep displacement curve and 
stability discriminant (in Chinese). J. Mountain Sci. 2000, 18, 440-444. 

10. Stiros, S.C.; Vichas, C; Skourtis, C. Landslide monitoring based on geodetically derived distance 
changes. J. Surv.Eng. 2004, 130, 156-162. 

11. Matsuura, S.; Asano, S.; Okamoto, T.; Takeuchi, Y. Characteristics of the displacement of a 
landslide with shallow sliding surface in a heavy snow district of Japan. Eng. Geol. 2003, 69, 15-35. 

12. Cornforth, D. Landslides in Practice: Investigations, Analysis, and Remedial/Preventive Options 
in Soils; John Wiley & Sons: Hoboken, NJ, USA, 2005. 

13. Atefi Monfared, K.; Rothenburg, L. Ground surface displacements and tilt monitoring for 
reconstruction of reservoir deformations. Int. J. RockMech. Mining Sci. 2011, 48, 1 1 13-1 122. 

14. Gili, J. A.; Corominas, J.; Rius, J. Using global positioning system techniques in landslide 
monitoring. Eng. Geol. 2000, 55, 167-192. 

15. Stramondo, S.; Saroli, M.; Tolomei, C; Moro, M.; Doumaz, F.; Pesci, A.; Loddo, F.; Baldi, P.; 
Boschi, E., Surface movements in Bologna (Po Plain — Italy) detected by multitemporal DInSAR. 
Remote Sens. Environ. 2007, 110, 304-316. 

16. Tomas, R.; Herrera, G.; Delgado, J.; Lopez-Sanchez, M.; Mallorqui, J.; Mulas, J. A ground 
subsidence study based on DInSAR data: Calibration of soil parameters and subsidence prediction 
in Murcia City (Spain). Eng. Geol. 2010, 777, 19-30. 

17. Dunnicliff, J. Measurement of deformation. In Geotechnical Instrumentation for Monitoring Field 
Performance', John Wiley & Sons: Hoboken, NJ, USA, 1993; pp. 199-249. 

18. Chen, Z. Landslides and engineered slopes. From the past to the future. In Proceedings of the 
Tenth International Symposium on Landslides and Engineered Slopes, Xi'an, China, 30 June-4 July 
2008; Volume 1. 

19. Bloyet, J.; Beghoul, N.; Ricard, Y.; Froidevaux, C. In situ test of a borehole extensometer. Rock 
Mech. Rock Eng. 1989, 22, 289-297. 

20. Bayoumi, A. On the evaluation of settlement measurements using borehole extensometers. 
Geotech. Geol. Eng. 2011, 29, 75-90. 

21. Corominas, J.; Moya, J.; Lloret, A.; Gili, A.; Angeli, M.G.; Pasuto, A.; Silvano, S. Measurement 
of landslide displacements using a wire extensometer. Eng. Geol. 2000, 55, 149-166. 



Sensors 2012, 12 



259 



22. Kalkani, E.C. Filtering probe inclinometer data to identify characteristics of slope movement. 
RockMech. Rock Eng. 1980, 13, 57-69. 

23. Simeoni, L.; Mongiovi, L. Inclinometer monitoring of the castelrotto landslide in Italy. J. Geotech. 
Geoenviron. Eng. 2007, 133, 653-666. 

24. Stark, T.; Choi, H. Slope inclinometers for landslides. Landslides 2008, 5, 339-350. 

25. Dowding, C.H.; Huang, F.C. Early detection of rock movement with time domain reflectometry. 
J. Geotech. Eng. 1994, 120, 1413-1427. 

26. Kane, W.; Beck, T.J. An alternative monitoring system for unstable slopes. Geotech. News 1996, 
14, 24-26. 

27. Chen, R.-P.; Chen, W.; Chen, Y.-M. TDR measurement system and the application of TDR in 
geoenvironmental engineering. In Advances in Environmental Geotechnics, Chen, Y., Zhan, L., 
Tang, X., Eds.; Springer: Berlin, Germany, 2010; pp. 153-162. 

28. Shentu, N.; Zhang, H.; Li, Q.; Zhou, H. Research on an electromagnetic induction-based deep 
displacement sensor. IEEE Sens. J. 2011, 77, 1504-1515. 

29. Cruden, D. A simple definition of a landslide. Bull. Eng. Geol. Environ. 1991, 43, 27-29. 

30. Cruden, D.; Varnes, D. Landslide types and processes. In Landslides — Investigation and 
Mitigation: Transportation Research Board Special Report 247; Turner, A.K., Schuster, R.L., Eds.; 
National Research Council: Washington, DC, USA, 1996; pp. 36-75. 

31. Highland, L.; Bobrowsky, P.T.; Survey, G. The Landslide Handbook: A Guide to Understanding 
Landslides; US Geological Survey: Reston, VA, USA, 2008. 

32. Nasar, S.A.; Xiong, G. Determination of the field of a permanent-magnet disk machine using the 
concept of magnetic charge. IEEE Trnas. Magn. 1988, 24, 2038-2044. 

33. Selvaggi, J.P.; Salon, S.; Kwon, O.M.; Chari, M.V.K. Calculating the external magnetic field 
from permanent magnets in permanent-magnet motors-an alternative method. IEEE Trans. Magn. 
2004, 40, 3278-3285. 

34. Selvaggi, J.P.; Salon, S.; Kwon, O.M.; Chari, M.V.K.; DeBortoli, M. Computation of the external 
magnetic field, near-field or far-field, from a circular cylindrical magnetic source using toroidal 
functions. IEEE Trans. Magn. 2007, 43, 1 153-1 156. 

35. Ravaud, R.; Lemarquand, G.; Lemarquand, V.; Depollier, C. Analytical calculation of the 
magnetic field created by permanent-magnet rings. IEEE Trans. Magn. 2008, 44, 1982-1989. 



© 2012 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article 
distributed under the terms and conditions of the Creative Commons Attribution license 
(http ://creati vecommons . org/licenses/by/ 3.0/). 



